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ABSTRACT 


Several  recent  papers  have  treated  the  high  efficiency  avalanche 
diode  using  a  simplified  analysis  which  successfully  explains  the  principal 
features  of  the  traveling  avalanche  zone,  formation  of  the  trapped  plasma, 
and  plasma  extraction.  These  analyses  have  assumed  square  wave  current 
excitation  so  that  all  phenomena  occur  under  constant  current  conditions. 

In  this  paper  it  is  shown  that  the  square  wave  current  excitation  requires 
negative  impedance  terminations  above  the  5th  harmonic.  Since  the  harmonic 
amplitudes  in  a  square  wave  decrease  only  as  l/n,  harmonics  as  high  as 
the  9th  and  11th  must  be  properly  terminated  to  make  the  analysis  realistic 
and  to  prevent  premature  avalanche.  We  have  analyzed  the  TRAPATT 
operation  using  the  same  physical  approximations  for  plasma  formation 
and  recovery  as  used  in  the  previous  square  wave  analyses.  However,  we 
have  assumed  a  trapezoidal  current  waveform  which  contains  primarily  first 
and  third  harmonics  with  negligibly  small  higher  harmonics.  In  this  analysis 
the  avalanche  zone  transit  and  portions  of  the  plasma  extraction  process 
occur  while  the  current  is  changing.  The  resulting  device  voltage  is  calcu¬ 
lated  as  a  continuous  analytical  function  of  time  except  during  the  plasma 
formation  interval  and  the  final  saturated  velocity  phase  of  the  plasma 
extraction  interval,  where  an  interpolated  voltage  waveform  is  obtained. 
Fourier  analysis  of  the  voltage  waveform  provides  the  device  efficiency 
as  well  as  the  necessary  circuit  termination  at  the  1st  and  3rd  harmonics. 
Realizability  conditions  on  the  diode  parameters  are  also  obtained. 
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AN  ANALYSIS  OF  FINITE  RISE  AND  FALL  TIME 
CURRENT  EXCITATION  OF  TRAPATT  DIODES 


I.  Introduction 

An  analysis  is  presented  for  the  operation  of  an  avalanche  diode 

oscillator  in  the  "trapped  plasma"  mode  when  excited  by  a  finite  rise-time 

12  3 

(trapezoidal)  current  waveform.  Previous  papers  *  *  analyzing  the  oper¬ 

ation  of  the  trapped  plasma  mode  have  detailed  the  traveling  ionization 
zone,  formation  of  the  trapped  plasma,  and  the  exit  of  the  plasma  under 
constant  current  driving  conditions.  Furthermore,  in  Refs.  1  and  2  approx 
imate  closed  form  analytical  solutions  have  been  presented  for  calculating 
a  complete  cycle  of  the  steady  state  voltage  waveform  assuming  an  ideal 
square  wave  current  drive.  The  analysis  using  a  square  wave  current 
provides  considerable  insight  into  the  details  of  the  trapped  plasma  cycle. 
This  paper  extends  the  analysis  to  handle  a  somewhat  more  realistic  cur¬ 
rent  waveform  with  finite  rise  and  fall  times.  The  design  information 
resulting  from  this  extended  analysis  provides  a  more  accurate  basis  for 
calculation  of  circuit  impedances  and  operating  frequencies  while  at  the 
same  time  relaxing  the  requirement  for  control  of  high  harmonic  termina¬ 
tions. 


II.  Rise  Time  and  Harmonic  Considerations  with  Square  Wave  Current 

Excitation 

Assume  the  analysis  of  Ref.  1  (hereafter  called  the  C1N  analysis 
after  Clorfeine,  Ikola  and  Napoli)  is  applied  to  the  diode  used  for  illustration 
in  that  reference.  The  oscillator  characteristics  are  given  in  Table  I.  The 
current  and  voltage  waveforms  are  shown  in  Fig.  1.  Fourier  analysis  of 
the  calculated  voltage  waveforms  results  in  diode  impedances  at  the  signi¬ 
ficant  (odd)  harmonics  as  shown  in  Fig.  2.  (The  current  waveform  has  no 
even  harmonic  content,  resulting  in  infinite  impedances  at  the  even  har¬ 
monic  frequencies). 
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TABLE  I 


Characteristics  of  TRAPATT  Oscillator  Using  CIN  Analysis  (Ref.  1) 


Diode  Type: 
n-layer  Doping: 
n-layer  Width: 

Saturated  Electron  Velocity: 

Saturated  Hole  Velocity: 

Low  Field  Electron  Mobility: 

Low  Field  Hole  Mobility: 

Area: 

Critical  Current  Density: 

Peak  Current  Density: 

Normalized  Current  Density: 

Pre -Avalanche  Electron  Concentration: 
Multiplication  Factor: 

Ionization  Constant: 

Dynamic  Breakdown  Field: 

Period: 

Frequency  of  Operation: 


Silicon  -  p  nn 

n  1 5  -3 

N ,  =  10  cm 

Cl 

W  =  8x10  ^cm 

7  -1 

v  =1.0x10  cm-sec 
ns 

7  -1 

v  =  1.0x10  cm-sec 
ps 

^  =  1220  cm^-v  1-sec  * 


*n 


UT 


2  - 1  - 1 
=  460  cm  -v  -sec 


-3  2 

A  =  3.  02  x  1 0  cm 

J  =  1.6x10^  A/cm^ 
c 

J  =  3.  52  x  10^  A/ cm^ 


j  =  j/jc  =  2.  2 

, n10  -3 

n  =10  cm 


o 

m  =  14.  5 

(y  =  3.  69  x  lO^cm  * 

E  =  3.  29  x  10^  v-cm 
c 

T  =  1 . 06  nsec 
f  =  942  MHz 


2 


T/2 


T 


t(sec) 


(a) 


TIME  (sec) 

(b) 

Fig.  1.  Waveforms  of  TRAP  ATT  oscillator  using  CIN  analysis 
and  diode  characteristics  of  Table  I.  (a)  Assumed  square  wave 
current  drive,  (b)  Calculated  voltage  response  with  linear  inter¬ 
polation  in  plasma  formation  interval. 


Fig.  2.  Diode  impedances  at  odd  harmonic  frequencies  for 
TRAP  ATT  oscillator  using  CIN  analysis  and  diode  charac¬ 
teristics  of  Table  I  (junction  area  A  =  3.  02  x  10"^cm^). 
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It  is  seen  in  Fig.  2  that  the  7th,  9th  and  11th  harmonics  do  not  represent 
negative  diode  resistances.  Thus,  if  the  diode  is  to  operate  very  nearly 
as  analyzed,  it  is  necessary  to  terminate  the  high  harmonics  in  negative 
resistances  to  provide  a  match.  In  other  words,  the  diode  must  be  pumped 
at  these  frequencies. 

In  Ref.  1,  this  problem  is  bypassed  by  neglecting  all  harmonics  above 
the  5th.  It  is  not  always  possible  to  do  this.  When  the  high  harmonics  are 
eliminated  from  the  square  current  waveform,  the  current  rise  and  fall 
times  increase.  When  this  happens,  the  trapped  plasma  cycle  may  be 
significantly  altered.  One  of  the  most  serious  consequences  of  increasing 
the  current  rise  and  fall  times  is  that  premature  ionization  can  occur 
either  because  the  current  does  not  shut  off  fast  enough  at  the  end  of  the 
plasma  extraction  period  or  because  it  does  not  turn  on  fast  enough  at  the 
beginning  of  the  charging  period. 

Premature  ionization  must  be  avoided  at  the  end  of  the  residual 
extraction  period  in  the  high  efficiency  cycle.  If  premature  breakdown 
occurs  at  the  end  of  plasma  extraction,  the  high-voltage  low-current 
portion  of  the  cycle  is  eliminated  and  the  efficiency  of  the  oscillator  is 
drastically  reduced.  Premature  ionization  during  the  charging  period  is 
not  as  serious  a  problem,  although  it  may  alter  the  plasma  field  and  plasma 
carrier  concentrations. 

The  limits  on  rise  and  fall  time  may  be  examined  quantitatively  by 

assuming  a  drive  current  with  finite  rise  and  fall  times.  Consider  the 

current  fall  time  first.  Assume  the  current  starts  to  fall  at  the  end  of 

plasma  extraction  (t  =  t^).  The  field  E  at  the  p  -  n  junction  at  this 

time  depends  on  the  peak  current  and  the  N^W  product.  These  are  chosen 

so  that  E  <  E  .  where  E,  is  the  static  breakdown  field.  E  continues  to 
w  b  b 

increase  during  the  current  fall  time  until  the  current  decreases  to  zero. 
The  final  field  at  the  junction  is  E^,  the  holding  field.  If  for  example,  the 
current  is  assumed  to  fall  linearly  according  to  the  equation 
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(1) 


J(t>  =  J  -  -  '»)•  *w  <  *  <  ‘w  +  t£ 

Then  the  holding  field  is  given  by: 

t  +  T, 

W  f  JT 

E  =  E  I1  f  J(t)dt  =  E  +  -=-=■  (2) 

h  w  e  w  w  2e  '  ' 


The  maximum  allowable  fall  time 
setting  Eh  =  Efe,  giving: 


t,  to  prevent  avalanche  is  obtained  by 

f  max  ^  1 


t  =  (E  -  E  )  ^ 

f  max  b  w  J 


(3) 


In  reality,  the  holding  field  must  be  less  than  E^  in  order  to  allow 

additional  time  for  the  current  to  rise  without  causing  low  current  ionization 

at  the  end  of  the  holding  period.  The  maximum  allowable  current  rise  time 

at  the  end  of  the  holding  period  T£max  calculated  in  the  same  way  as 

t  assuming  the  current  rises  linearly  from  zero  to  J  in  time  T  . 

f  max  &  /  r 


T 

r  max 


<Eb  -  Eh>x 


(4) 


Equal  time  intervals  may  be  provided  for  current  turn-on  and  turn¬ 
off  by  choosing 


E 


h 


E 

w 


+  Eb 


2 


(5) 


Then  the  single  critical  time  constant  holds  for  rise  and  fall 

t  =  (E  -  E  J  7  (6) 

max  '  b  w  J 

The  occurrence  of  low  current  ionization  during  the  charge  period 
does  not  alter  the  high  efficiency  cycle  in  a  first  order  way  unless  it  occurs 
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before  the  current  has  reached  the  critical  level  J  for  formation  of  a 

c 

traveling  ionization  zone,  where  =  qN^v^.  If  the  current  reaches  the 
level  before  the  junction  field  charges  to  E^,  the  analysis  can  be  carried 
out  as  in  the  ideal  square  wave  case  with  appropriate  modifications  to  account 
for  the  changing  current  during  the  transit  of  the  ionization  wave. 


max 

amp/ 


A  value  can  be  calculated  for  the  maximum  allowable  transition  time 

3 

for  the  sample  diode  of  Table  I.  For  this  diode,  using  J  =  3.  52  x  10 


cm^,  E  =  3.  29  x  10^  v/  cm,  and  E  =  2. 85x10^  v/  cm: 
b  w  ' 


t  =  1.3x10  ^sec  (7) 

(The  value  of  E^  used  here  is  the  dynamic  breakdown  field  E^  of  Ref.  1.  ) 

This  time  may  now  be  compared  with  the  rise  time  obtained  by 
deleting  all  the  high  order  Fourier  components  in  the  analysis  of  an  ideal 
square  wave  with  the  fundamental  frequency  given  in  Table  I.  Figure  3 
illustrates  the  normalized  leading  edge  waveforms  for  Fourier  synthesized 
square  waves  with  partial  sums  varying  from  1st  harmonic  alone  to  1st 
through  11th  harmonics.  The  normalized  times  required  for  the  various 
waves  to  increase  from  -1  to  +1  are  tabulated  in  Table  II.  Also  tabulated 


are  the  actual  rise  times  if  the  period  is  1. 06  nsec  as  indicated  in  Table  I. 

TABLE  II 

Rise  Times  for  Waveforms  of  Fig.  3 


Highest  Harmonic 

Rise  Time  Normalized 

Actual  Rise  Times  for 

Content 

to  Period 

Period  T  =  1.06  nsec 

(sec) 

i 

0.  288 

3.  04  x  10’10 

3 

0.  151 

1 .  63  x  10  ^ 

5 

0.  101 

1. 07  x  10’10 

7 

0.  077 

8.12x10  *  * 

9 

0.  061 

6.49  x  10’U 

1 1 

0.  051 

5.  42  x  lO-11 
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NORMALIZED  AMPLITUDE 


Fig.  3.  Leading  edge  waveforms  for  partial  Fourier  sums  of  an 
ideal  square  wave  with  normalized  peak-to-peak  amplitude  of  2. 
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A  comparison  of  this  tabulation  with  Tmax  indicates  that,  despite  the 
inclusion  of  1 1  harmonics,  the  maximum  allowable  current  rise  time  is 
exceeded  by  about  a  factor  of  4. 

If  only  the  5th  harmonic  is  included,  t  is  exceeded  by  almost 
7  max  7 

a  factor  of  10.  t  may  be  increased  without  1st  order  changes  in  the 

high  efficiency  cycle  or  analysis.  If  the  N^W  product  is  decreased,  the 

field  at  the  end  of  plasma  extraction  is  decreased  and,  according  to 

Eq.  (6),  t  is  increased.  Decreasing  J  also  increases  t  .  How- 

ever,  when  using  realistic  diode  parameters  and  current  densities  the 

square  wave  analysis  does  not  generally  result  in  j  sufficiently  large 

max 

to  permit  the  elimination  of  harmonics  higher  than  the  5th  in  the  driving 
waveform.  Finite  rise  and  fall  times  are  required  for  reduced  harmonic 
content  in  the  current  waveform.  Thus  the  square  wave  analysis  must  be 
extended  to  handle  plasma  formation  and  plasma  extraction  under  variable 
current  conditions. 

III.  Trapezoidal  Waveform 

It  is  possible  to  control  the  onset  of  avalanche  and  simulta¬ 
neously  reduce  the  harmonic  content  of  the  current  waveform  by  increasing 
the  rise  and  fall  times  of  the  current  waveform  and  allowing  plasma  forma¬ 
tion  to  occur  while  the  current  is  rising  and  plasma  extraction  to  continue 
while  the  current  is  falling.  For  the  purpose  of  analysis,  the  square  wave 
current  is  replaced  with  the  ideal  trapezoidal  waveform  of  Fig.  4  (solid 
line).  Several  trapezoidal  waveforms  have  been  investigated.  The  wave¬ 
form  shown  in  Fig.  4  has  a  peak-to-peak  rise  time  which  is  exactly  T/5, 
where  T  is  the  period.  This  waveform  has  the  special  property  that  its 
5th  harmonic  is  identically  zero. 
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NORMALIZED  AMPLITUDE 


Fig.  4.  Trapezoidal  driving  current  waveform  with  partial  Fourier  sums. 


The  amplitudes  of  the  harmonics  of  this  waveform  decrease  approximately 
2 

as  l/n  .  The  trapezoidal  waveform  contains  only  odd  harmonics  (a  conse- 

T 

quence  of  the  fact  that  J(t  +  — )  =  ” J( t) )  and  when  the  origin  is  properly 
chosen  the  trapezoidal  waveform  can  be  either  an  even  or  odd  function  of 
time.  Fourier  analysis  of  this  waveform  gives 

20J 

J(t)  =  — 2^(0.  sin  COt  +  0.  1056  sin3uot  -  0.  0194  sin7u)t-.  .  . )  (8) 

TT 

It  is  seen  that  the  amplitude  of  the  7th  harmonic  of  this  waveform  is  only 
3.  3%  of  the  fundamental,  so  that  the  elimination  of  harmonics  above  the 
third  in  this  waveform  has  little  effect  on  its  shape.  This  can  be  seen  by 
examining  the  dotted  and  opened  curves  of  Fig.  4  which  are  plots  of  partial 
Fourier  sums  of  this  waveform  including  respectively  harmonics  through 
the  third  and  seventh.  It  is  clear  that  the  third  harmonic  partial  sum  is  a 
good  approximation  to  the  ideal  trapezoidal  waveform  and  that  the  rise  time 
is  changed  only  slightly  as  the  seventh  harmonic  is  added. 

IV.  Diode  Response  to  Trapezoidal  Current  Drive 

The  solution  for  the  diode's  response  to  this  trapezoidal  current 
excitation  is  based  on  the  same  approximations  for  the  ionization  coefficient 
OL  that  were  used  in  Ref.  1.  The  analysis  proceeds  along  somewhat  different 
lines  in  that  the  one  cycle  of  the  current  waveform  is  divided  into  four  seg¬ 
ments  (rise,  hold  at  maximum  current,  fall,  hold  at  zero  current)  as  opposed 
to  the  two  segments  of  the  square  wave  current  drive.  The  resulting  voltage 
waveform  still  retains  the  same  basic  features  as  that  arising  from  a  square 
wave  drive  except  for  a  certain  amount  of  smoothing  corresponding  to  the 
smoothing  of  the  current  waveform.  The  trapezoidal  current  waveform  is 
divided  into  ten  equal  time  increments,  At.  The  rise  and  fall  times  have 
a  duration  of  2At  and  the  flat  portion  of  the  waveform  has  a  duration  of 
3At.  There  is  some  latitude  available  in  specifying  at  what  point  in  the 
current  cycle  each  of  the  critical  transitions  in  the  voltage  waveform  is 
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assumed  to  take  place.  The  operating  cycle  to  be  outlined  below  is  one 
of  several  that  could  be  chosen.  The  various  transitions  in  the  cycle  are 
chosen  either  (a)  for  the  sake  of  analytical  convenience,  (b)  to  provide 
the  necessary  phase  relationship  between  voltage  and  current  to  assure 
that  all  significant  harmonics  can  be  terminated  in  positive  current,  or 
(c)  to  assure  that  the  field  at  the  p-n  junction  remains  below  the  static 
breakdown  field  at  all  times.  This  analysis  has  been  developed  with 
sufficient  generality  to  allow  the  rise  time  of  the  trapezoidal  current  wave¬ 
form  and  starting  time  of  residual  extraction  to  be  varied.  These  para¬ 
meters  have  been  selected  in  the  present  discussion  to  provide  proper 
phasing  and  to  prevent  premature  avalanche  with  a  particular  example 
treated  here.  Any  operating  cycle  is  valid  which  obeys  the  diode  physics 
and  which  leads  to  self-consistent  results  and  a  realizable  set  of  operating 
impedances.  The  physical  procedure  for  changing  from  one  realizable 
operating  cycle  (or  mode)  to  another  is  then  to  alter  the  current  waveform 
by  retuning  the  external  oscillator  circuitry  and/or  by  changing  the  current 
amplitude. 

In  the  process  of  searching  for  realizable  current  waveforms,  it 
became  apparent  that  the  trapped  plasma  oscillator  is  not  capable  of 
functioning  over  a  wide  range  of  current  rise  times  or  amplitudes.  The 
amplitude  and  rise  time  must  be  constrained  to  lie  within  relatively  narrow 
values  to  insure  proper  phasing  to  prevent  unwanted  breakdown  transients. 
This  observation  is,  of  course,  totally  consistent  with  experimental  results. 
One  of  the  consequences  of  specifying  a  ratio  of  current  rise-time  to  period 
and  a  specific  time  for  start  of  avalanche  is  that  the  depletion  layer  width, 
doping  concentration,  and  current  density  can  no  longer  be  specified  inde¬ 
pendently.  For  analytical  convenience,  the  average  current  density  (J  ) 

a 

during  plasma  formation  and  the  base  doping  density  are  chosen  as 
independent  variables,  and  the  width  and  operating  frequency  become  depen¬ 
dent  variables.  The  analysis  assumes  an  ideal  p+ n  n+  diode  with  uniformly 
doped  n  region.  In  the  following  description  the  terminology  used  to  describe 
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the  various  parts  of  the  cycle  will  be  consistent  with  that  used  in  Ref.  1. 

The  analysis  determines  values  of  electric  field  and  voltage  as  continuous 
functions  of  time  during  most  of  this  cycle  except  for  plasma  formation 
and  residual  plasma  extraction  during  which  a  continuous  time  dependent 
voltage  waveform  is  generated  by  interpolation.  In  addition,  values  of 
base  width  W  and  period  T  are  also  obtained  in  the  analysis.  This 
description  of  the  analysis  will  be  divided  into  sections  corresponding  to 
the  lettered  regions  shown  in  Fig.  5.  The  field  profiles  corresponding  to 
these  regions  are  plotted  in  Fig.  6. 

A.  Charging  Period:  Region  A  (0  £  t  £  At) 

The  cycle  begins  at  a  time  when  the  maximum  field  in  the  depletion 
region  is  less  than  the  static  breakdown  field  but  large  enough  to  cause 
punch  through  into  the  n*  region.  At  this  point  the  diode  voltage  is  defined 
as  V  and  the  peak  field  is  E  ,  both  of  which  are  to  be  determined.  In 
the  absence  of  ionization  there  are  no  free  carriers  in  the  depletion  region 
so  only  displacement  current  flows: 

J(t>  =  «§f  (9) 


It  can  be  seen  from  Fig.  6  that  during  this  interval  J(t)  is  given  by 


(10) 


where  Jq  is  the  maximum  value  of  the  trapezoidal  current  waveform  and  is 
yet  to  be  determined.  At  is  specified  such  that  the  electric  field  increases 
from  to  a  value  E  just  as  the  current  reaches  the  level  J q/2,  giving 


At 


-  E  ) 
m 


(ID 
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(0) 


(b) 


Fig.  5.  Waveforms  of  TRAP  ATT  oscillator  using  trapezoidal  drive, 
(a)  Assumed  drive  current,  (b)  Calculated  voltage  response  with 
linear  interpolation  in  plasma  formation  interval. 
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©  © 


Fig.  6.  Approximate  field  profiles  in  n-layer  of  TRAPATT 
oscillator  using  trapezoidal  current  drive.  Each  diagram 
illustrates  corresponding  lettered  interval  in  previous  figure. 
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is  the  peak  field  attained  at  the  p+n  junction.  Both  E  ^  and  Jq  can 
be  calculated  once  J  ,  the  average  current  density  during  plasma  formation, 

3L 

has  been  specified.  will  be  discussed  in  detail  in  the  next  section. 

The  time  dependent  voltage  during  this  period  is: 

t  2 

V(t)  =  Vm  +  Vpk(^),  (0<ts  At)  (12) 


where 


V 


pk 


WE 


Vk 


qN  ,W 
d 


(13) 


B.  Plasma  Formation  Period:  Region  B  (At  £t  <At  +  t^) 

Plasma  formation  occurs  if  the  current  density  exceeds  a  critical 
value,  J  =  qN^v  before  the  junction  field  has  reached  its  dynamic  break¬ 
down  level  E^.  When  this  happens  an  ionization  zone  sweeps  across  the 
depletion  layer  with  velocity  v^  leaving  in  its  wake  a  dense  plasma  of  holes 
and  electrons  with  concentrations  n^  and  p^.  Expressions  for  the  resulting 
particle  concentrations  and  fields  are  given  in  the  literature  based  on  the 
assumption  of  a  constant  zone  velocity  which  implies  a  constant  diode  cur¬ 
rent  J(t)  during  the  time  t^  that  the  zone  is  in  the  depletion  layer.  In  this 
analysis, the  current  changes  during  zone  transit.  However,  it  is  valid  to 
use  the  constant  current  approximation  if  the  zone  transit  occurs  in  a  time 
interval  which  is  small  compared  to  the  rise -time  of  the  current  excitation. 
The  zone  transit  time  is 

w  +  6. 

t£  =  — r- £qNd  <14> 

a 


where  6  is  the  plasma  sheath  width  and  J  is  the  average  current  density 
s  a 

during  zone  transit.  The  quality  of  the  approximation  depends  on  the  ratio 
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(15) 


In  most  cases  is  approximately  0.  5,  indicating  that  the  zone  tran¬ 
sit  occupies  one-fourth  of  the  rise-time  and  the  approximation  is  not  unrea 
sonable. 


In  order  to  make  use  of  the  equations  for  plasma  formation  previously 

presented  in  the  literature  it  is  necessary  to  know  the  current  density  during 

plasma  formation.  Thus,  in  this  analysis  J  assumes  the  role  of  an  input 

a 

parameter  in  the  analysis  along  with  N^.  Given  and  we  may  proceed 
directly  to  calculate  all  of  the  plasma  parameters.  We  will  follow  the  anal¬ 
ysis  presented  in  Ref.  1  and  briefly  review  here  the  sequence  of  expressions 
derived  there  for  computing  the  plasma  parameters. 

First,  we  define  a  normalized  average  current  density 


Ja  qN.v 


d  ns 


(16) 


Then  the  peak  field  at  the  junction  is  given  by 

qN 

E  =E  + - -(j  -l)jen(Y-l)  (17) 

pk  c  a  e  a  ' 


where 


Y  = 


Ja+r\Nd 


1  +  r  /  n 


(18) 


O'  =  2.4x10  ^  exp(  -1.6  x  1 0^/E  cm  * 

(the  ionization  coefficient  for  silicon) 
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r  =  v  /v  (the  ratio  of  saturated  hole  and  electron  (19) 

ps  ns  velocities) 


-1 


E  = 
c 


InZ  +  1 


(the  dynamic  breakdown  field  (20) 
U.  6  x  10U  Ep^  f  at  the  junction) 


10  -3 

n  =10  cm  (the  electron  concentration  in  the  depletion 
layer  before  ionization) 


The  plasma  field  is 

E  = 


(j  )2v 
XJa'  ns 


P  =  < ja  +  r)(Mn  +  /*p)(Ps/Nd)  +  ia<A*n  -  Mp) 


(21) 


where 


2  - 1  - 1 

fji  =  1220  cm  volt  sec  (the  electron  mobility) 


2  "1  “1 

fJL p  =  460  cm  volt  sec  (the  hole  mobility) 


ps  =  ^  y1  ^  —  J  n^e111  (the  sheath  hole  density) 


(22) 


m  =  In  (y  in  Y)  (the  multiplication  factor) 


(23) 


The  sheath  width  is 


eE 


6_  = 


s  q(n  -  p  -  n.) 

s  *s  d' 


(24) 


where 


ng  =  nQem  (the  sheath  electron  density) 


(25) 
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The  plasma  electron  and  hole  densities  are  respectively 


n 

P 


j  v 

Ja  ns 


l)v  n 
ns 


-  U  E 
"n 


s 

P 


PP  = 


M  v  +  v 
Ja  ns  ps 

j  v  +  ii  E 
Ja  ns  rp  p/ 


(26) 


(27) 


Although  all  of  the  plasma  parameters  can  be  calculated  at  this  point,  the 
calculation  of  depletion  layer  width  W  and  zone  transit  time  t^,  which 
depends  on  W  must  be  deferred  to  a  later  section. 

The  time  dependence  of  voltage  during  plasma  formation  is  relatively 
complex.  We  follow  the  procedure  of  Ref.  1  and  approximate  it  with  a 
linear  decrease  in  voltage  during  the  zone  transit.  This  approximation  is 
reasonable  since  the  zone  transit  time  t^  is  a  sufficiently  small  fraction 
of  the  total  TRAPATT  period  that  a  more  detailed  analysis  would  have 
negligible  affect  on  the  efficiency  and  impedance  calculations.  Neglecting 
the  effects  of  the  drift  of  the  low  velocity  electrons  during  the  ionization 
zone  transit  and  the  transient  effects  as  the  zone  exits  the  depletion  layer 
a  uniform  plasma  fills  the  depletion  layer  after  the  zone  transit  is  complete. 
The  resulting  minimum  voltage  is  then 


V  =  WE  (28) 

P  P 


Hence  the  time  dependent  voltage  during  zone  transit  is 

V(t)  =  VPk  +  (VP  ’  vPk}  (^  *  t  <;  At  +  tf)  (29) 


C.  Rising  Current  Plasma  Region:  Region  C  (At  +  t  s  2At) 

After  the  plasma  has  formed,  the  current  density  continues  increas¬ 
ing  until  it  reaches  Jq.  The  current  and  electric  field  are  ohmically  related 
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in  the  plasma-filled  region  so  that  the  electric  field  also  increases  linearly 

during  this  time.  The  peak  current  density  Jq  can  be  shown  to  be  related 

to  J  as 
a 


4J 

J  =  - — 

Jo  bf  +  2 


(30) 


and  the  electric  field  in  the  plasma  at  the  end  of  this  period  is 


E 

e 


4E 

_ E_ 

b£  +  2 


(31) 


An  exact  calculation  of  the  voltage  at  the  end  of  this  period  is  difficult. 
However  the  voltage  is  negligible  compared  to  the  breakdown  voltage  so 
that  simple  approximations  may  be  used  with  little  effect  on  the  final 
results.  Thus,  during  this  interval  it  is  again  assumed  for  purposes  of 
calculating  voltage  that  no  plasma  drift  occurs  and  the  depletion  layer 
remains  completely  filled  with  plasma.  The  resulting  end  voltage  is  then 

Ve  =  WEe  (32) 

and  the  time  dependent  voltage  expression  is 

V(t)  =  vp  +  <Ve  -  vp)  (~tA.*  tf 2At)  (33) 

In  order  to  determine  W  it  will  be  necessary  to  know  how  far  the 
holes  or  electrons  move  during  the  entire  plasma  extraction  cycle.  The 
distance  moved  by  the  holes  during  this  portion  of  the  cycle  is 

2At  At  -  t 

d  =  J  v(t)dt  =  J*  V  E'(t)dt  (34) 

Pl  At  +  t  0  P 


20 


where 


E'(t)  =  E  +  (E 
p  e 


giving 


d 

Pi 


P„ED  l>  -  bf)(6  +  b£) 
2  ^  2  +  b 


(35) 


(36) 


The  distance  moved  by  the  electrons  in  the  same  time  is 


d 

ni 


Pi 


(37) 


D.  Constant  Current  Plasma  Region:  Region  D  (2At  ^  t  ^  5At) 


As  holes  drift  out  of  the  depletion  region  they  leave  behind  a 
region  of  net  negative  charge  near  the  n^n  boundary  as  shown  in  Fig.  7a. 
A  similar  region  of  net  positive  charge  develops  near  the  p^n  boundary. 
The  departure  from  neutrality  in  these  two  regions  results  in  the  approx¬ 
imate  field  profile  shown.  We  follow  Ref.  1  in  approximating  the  field 
and  concentration  transitions  as  abrupt.  This  is  illustrated  in  Fig.  7b. 

The  velocity  in  the  non-neutral  regions  is  assumed  to  be  saturated.  The 
following  expressions  for  the  ’’residual’1  particle  densities  in  the  saturated 
velocity  regions  result  from  the  principle  of  current  continuity: 


n 

r 


(U  l|i  )E  n 
^n  ^p7  e  p 

v  +  JU  E 
ns  e 


Pr 


(H  +  n  )E  p 
n  p  erp 

v  +  u  E 
ps  'n  e 


(38) 

(39) 
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Fig.  7.  Electric  field  and  free  carrier  concentration 
profiles  during  plasma  extraction  period,  (a)  Actual 
profiles,  (b)  As  sumed  profiles. 
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To  calculate  the  voltage  during  this  portion  of  the  cycle  we  must 
integrate  the  changing  electric  field  profile  of  Fig.  7b.  The  result  is 
the  following  expression 

V(t)  .  EeW  +  1/2  jsn  [dpi  +  dp2(^)]2+  Sp[dni  +  dn2(^)]2j 

(40) 

where 

S  =  a(n  -  N,),  S  =  a(N,  +  p  )  (41) 

n  e  '  n  d7  p  e  d  *r'  '  ' 

are  the  field  gradients  in  the  residual  regions,  and 

d  _  =  3u  E  At,  d  _  =  3u  E  At  (42) 

p2  Pp  e  ’  n2  *n  e  '  ' 

are  the  total  distances  moved  by  the  trapped  carriers  during  the  constant 
current  phase. 

E.  Falling  Current  Plasma  Region:  Region  E  (5At  £  t  5At  +  t  ) 

a 

In  order  to  prevent  premature  avalanche  at  the  end  of  the 

plasma  extraction  process  it  is  necessary  to  allow  the  current  to  begin  to 

fall  before  the  trapped  plasma  disappears.  The  beginning  of  the  final 

saturated  velocity  plasma  recovery  phase  is  assumed  to  occur  t  seconds 

after  the  current  starts  to  fall.  For  a  given  set  of  diode  parameters  and 

drive  current  level  there  is  only  a  limited  range  of  realizability  for  t  . 

a 

If  t  is  too  short,  the  n-region  is  depleted  of  carriers  early,  allowing 
a 

the  junction  field  to  charge  to  breakdown  before  the  current  falls  to  zero. 

If  t  is  long,  the  free  carriers  may  not  be  cleared  out  of  the  n-layer  before 
a 

current  falls  to  zero.  In  this  analysis  it  was  found  that  specifying  t  =  0.  4At 

a 

was  sufficient  to  prevent  premature  breakdown  and  to  allow  time  for  all  the 
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residual  holes  and  electrons  to  drift  out  of  the  n-layer  before  the  current 
reaches  zero. 

The  physics  of  this  part  of  the  plasma  extraction  period  is  compli¬ 
cated  by  the  falling  current.  The  variation  of  the  carriers  and  the  electric 
field  in  this  interval  are  shown  in  Fig.  8.  In  this  figure  the  magnitude  of 
the  plasma  field  is  expanded  out  of  proportion  to  the  field  in  the  saturated 
velocity  regions.  Fig.  8a  shows  the  field  and  carrier  profiles  just  as  the 
current  starts  to  drop.  There  is  still  a  fraction  of  the  n-layer  filled  with 
the  plasma  concentration.  As  the  current  drops  the  field  in  the  ohmic 
plasma  region  drops  also.  As  a  result  the  residual  hole  and  electron 
densities  in  the  saturated  regions  also  decrease  with  time.  This  is  shown 
in  Fig.  8b,  where  it  is  clear  that  the  most  recently  generated  residual 
carriers  (nearest  the  plasma  region)  have  a  lower  concentration  than  the 
older  residual  carriers.  The  saturated  carrier  velocity  is  about  5  times 
the  plasma  velocity  so  that  the  reduced  carrier  concentrations  quickly 
spread  out  toward  the  boundaries  of  the  n-layer.  In  Fig.  8b  the  holes 
produced  as  the  current  first  started  to  fall  have  traversed  about  half 
the  distance  to  the  pfn  junction.  As  a  result  of  this  non-uniform  carrier 
concentration,  the  electric  field  profile  becomes  slightly  nonlinear  as 
shown.  In  Fig.  8c  this  process  has  progressed  far  enough  so  that  the 
entire  region  of  saturated  velocity  holes  is  non-uniformly  distributed. 

Fig.  8d  shows  the  field  and  carrier  profiles  at  the  instant  when  the  trailing 
edge  of  the  hole  packet  coincides  with  the  trailing  edge  of  the  electron 
packet  and  the  plasma  has  just  disappeared.  This  is  the  beginning  of  the 
saturated  velocity  "residual"  extraction  period. 

We  will  now  proceed  to  calculate  the  carrier  and  field  profiles  as 
a  function  of  time  during  this  interval  throughout  the  portion  of  the  n-layer 
which  is  filled  with  residual  holes.  The  electron  filled  portion  may  be 
analyzed  in  the  identical  manner  with  appropriate  variable  substitutions. 
The  end  result  is  an  expression  for  voltage  as  a  function  of  time  in  this 
interval.  During  this  period  the  field  in  the  (shrinking)  plasma  region 
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18-4-12955 


Fig.  8.  Field  and  carrier  profiles  during  falling  current  portion 
of  plasma  recovery  cycle  from  (a)  start  of  current  fall  to  (d)  end 
of  plasma  extraction. 
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is  given  by 


E  (t)  =  E  (1  - 

p  e 


(43) 


The  distance  moved  (to  the  left)  by  the  right  hand  edge  of  the  plasma  region 
is  the  distance  moved  by  the  trailing  edge  of  the  electron  packet. 


yn  = 


t  t 

y  y 

v  dt  =  u,  P  E  (t)dt 
J  np  'nJ  p'  ' 


jU  E  2 

U  E  t  - 
n  e  y  4At  y 


(44) 


solving  for  t 


t  (y  ) 
y'  7n' 


t  =  2  At 

y 


n 


u.  E.At 
^n  f 


(45) 


and  substituting  into  E^(t)  to  obtain  E^( yn) 


Ep<y„>  =  Ee'/1  -  ft  E,At 
r  'n  i 


n 


(46) 


This  is  the  field  which  produces  the  initial  residual  hole  concentration 

p  (y  )  at  the  moving  right  hand  edge  of  the  plasma  region*  p  (y  )  is 

o  o 

related  to  E  (y  )  by 
p  wn'  J 


5  (y  ) 

r  wn' 
o 


p  (/i  +  fJL  )E  (y  ) 

p'n  'p7  pwn/ 

v  +  u  E  (y  ) 
ps  pwn 


(47) 


This  relation  may  be  linearized  by  expanding 


E  (y  )  =  E  +  AE  (y  ) 
pwn'  e  pwn' 


(48) 
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where  E  is  the  plasma  field  at  the  start  of  the  current  fall,  and  AE  (y  ) 
e  r  pwn 

is  the  change  in  E  (y  ).  Since  AE  (y  )  «v  we  obtain 
&  p  7n  pwn'  ps 


p  (U  +  U  ) 

Pr  (vn)  -  JP  ».(y.>  -  c„Ejyj 


v  +  n  E  p'  'n' 
ps  *n  e  r 


P  P  n 


(49) 


The  residual  holes  drift  to  the  right  with  saturated  velocity  v^g.  At  anY 
time  the  distribution  of  holes  to  the  right  of  the  righthand  plasma  edge  is 


y  (t) 

p  (z  ,  t)  =  C  E  (y  (t) - 77TT - 777  z  ) 

*r'  n  '  p  pwn'  '  yn(t)  +  v„c(t)  n 


(50) 


ps' 


where  is  referenced  to  a  moving  coordinate  system  whose  origin  is  the 
righthand  edge  of  the  plasma  zone,  increases  positively  to  the  right  as 

shown  in  Fig.  8c.  The  field  in  this  region  of  nonuniformly  distributed 
holes  may  be  obtained  by  integrating  Poisson's  equation 

z 


n 


E  (z  ,  t)  =  H  (*  [N,  +  p  (z'  ) ] dz ' 

p'  n’  '  e  n  d  *V  n'  n 


(51) 


giving 


E  (z  ,t)  =  a 
pv  n  7  e 


r  3/2  3/2  ) 

NdZn  +  CpEe  3b  [(a(t)  +  b -  aW  ][  {52) 


where 


a(t)  =  1  -  TT  +  * 


At  4At^ 


(53) 


and 


b(t)  = 


At  "  4At2 


U  E.At  (-£-  -  -rrr?)  +  v  t 
pn  f  'At  4At^'  ps 


(54) 
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The  voltage  drop  across  the  region  is  then  obtained  by  integrating 


d 

n 

V  (t)  =  f  E  (z  ,  t)d  z 
P  J  0  p  n 


n 


N,d 
d  n 


3/2 


-  C  E 


2a(t) 


+  C  E 


C(a(t)  +  b(t)dn) 


p"e  3b(t)  '  '"p  ei5b2(t) 

5/2  5/2.  ) 

~  a(t)  ]  f 


(55) 


where  d  ,  the  width  of  the  non-uniform  distributed  residual  hole  region, 
n 

is  measured  from  the  edge  of  the  plasma  (z^  =  0). 

If  the  edge  of  this  region  has  not  drifted  to  the  p-n  junction,  i.  e.  , 

t<  d  /v  ,  (where  d  is  the  total  width  of  the  residual  hole  saturated 
nc  ps  nc 

velocity  region  at  the  start  of  the  current  fall),  then  the  voltage  in  the 
saturated  region  includes  another  term  Vx  which  is  the  integral  of  the 
linear  field  in  the  region  adjacent  to  the  junction. 

V  =  E  (d  (t),t)  .  (d  -  v  t)  +  0.  5S  (d  -  v  t)2  (56) 

x  p'  n'  "  '  '  nc  ps  '  p'  nc  ps  '  '  ' 


where 


dn<‘>  * 


yn(l> 


V  t 

ps 


(57) 


Expressions  analogous  to  (55)  and  (56)  can  be  obtained  for  the  saturated 
velocity  residual  electron  region.  The  only  remaining  component  of  voltage 
is  then  given  by  the  integral  of  the  initial  plasma  field  over  the  entire  n- 
layer  width  W  E^fl  -  2At^‘ 

During  this  time  interval  it  is  necessary  to  keep  track  of  the  peak 
field  reached  at  the  junction  to  assure  that  premature  breakdown  does  not 
occur.  The  rate  of  change  of  field  at  the  junction  depends  both  on  the 
external  applied  current  and  the  drift  current  of  residual  holes. 
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where  p  .(t)  is  the  instantaneous  hole  concentration  at  the  junction.  From 
^  3 

this  expression  it  can  be  seen  that  the  junction  field  will  reach  a  maximum 
and  start  decreasing  when  J(t)  has  decreased  sufficiently  to  make  the  right" 
hand  side  of  Eq.  (58)  equal  to  zero;  i.  e.  ,  when 

Jo<‘  -2St»  =  <59> 

(assuming  this  occurs  before  the  residual  carriers  are  completely  cleared 
out  of  the  n-layer). 

If  this  maximum  is  reached  before  Prj(t)  begins  to  decrease,  the 
time  of  occurrence  and  the  magnitude  of  the  maximum  junction  field  are 
easily  calculated. 

The  time  of  occurrence  measured  from  start  of  current  fall  is 

J 

t.  =  2At  [1  -  y-  ]  (60) 

J  o 


where 


J 


r 


qprv 


ps 


(61) 


and  p^  is  the  initial  residual  hole  concentration  obtained  during  the  con¬ 
stant  current  phase. 

The  maximum  field  at  the  junction  is 


J  At 
o 


Jr\2 


E.  =  E-  + 
imax  fp 


e 


il 


(62) 


where  E.  is  the  junction  field  at  the  start  of  the  current  fall  time.  E. 

fp  J  imax 

is  a  relative  maximum  only,  since  the  field  starts  rising  again  when  the 
residual  extraction  period  is  complete. 

If  the  maximum  field  occurs  after  the  residual  hole  concentration 
at  the  junction  starts  to  decrease,  the  maximum  field  computation  is  consid¬ 
erably  more  complicated.  This  computation  will  not  be  included  here  be¬ 
cause  it  is  not  needed  in  the  example  to  be  discussed  below. 

F.  Residual  Extraction  Region:  Region  F  (5At  +  t  ^t^5At+t  +t  ) 

cL  3.  IT 

When  the  trailing  edges  of  the  unsaturated  plasma  and  hole 
concentrations  coincide  as  shown  in  Fig.  8d,  the  residual  plasma  extraction 
period  begins.  The  field  now  exceeds  the  saturation  level  everywhere  in 
the  n-layer.  In  silicon  with  equal  saturation  velocities  for  holes  and  electrons 
the  remaining  time  required  to  completely  sweep  the  depletion  region  clear 
of  carriers  is  the  time  required  to  remove  the  residual  holes  (which  have 
the  greater  distance  to  travel. ) 
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I U  W 
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(63) 
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Exact  calculation  of  a  time  dependent  voltage  during  this  time  interval  is 
difficult.  For  this  analysis  it  suffices  to  interpolate  between  the  beginning 
and  end  voltages. 

The  field  at  the  end  of  residual  extraction  is  obtained  by  startigg 

from  the  minimum  field  at  the  beginning  of  residual  extraction  E  (1 

e  a 

and  following  it  as  it  charges  to  a  value 


E  ' 
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where 


b 


r 


(66) 


at  the  end  of  the  residual  extraction  time  period. 

From  Fig.  6f  it  is  seen  that  the  junction  field  at  this  time  is  then 

qN,  M  W 
d 


E  =  E  1  + 
n  n  e  11  +  U 


(67) 


The  voltage  at  the  end  of  residual  extraction  is  simply 

rZ 


V  =  E  W 
n  n 


qNdW 
2e 


(68) 


G.  Final  Charging  Region:  Region  G  (5 At  +  t  +  t  s  t  s  7 At) 

a  r 

In  this  time  interval  all  carriers  have  been  cleared  out  of 
the  n-layer  and  the  field  continues  charging  under  the  influence  of  the 
decreasing  current.  We  are  now  able  to  calculate  E  ,  the  junction  field 
at  the  end  of  the  current  fall  time,  in  the  same  way  E^  was  calculated, 
giving 


E  = 
m 


E 

e 


qN ,  a  W 
2 _ d_ 

e  u  +  u 

^n  rp 


(69) 


The  time  dependent  voltage  during  this  period  is  given  most  simply 
as  a  function  of  E 

n 
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(70) 


V(t)  =  W<E  + 
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H.  Zero  Current  Holding  Region:  Region  H  (7 At  ^  t  ^  lOAt) 

In  this  time  interval  the  n-layer  is  depleted  of  carriers  and 
no  current  flows.  Hence  the  field  profile  remains  unchanged  from  its 
value  at  the  end  of  interval  G.  At  the  end  of  this  interval  the  TRAPATT 
cycle  is  complete. 

I.  Determination  of  t^  and  W 

At  this  point  in  the  analysis  all  of  the  important  timing  and 
electric  field  relations  have  been  derived.  However,  these  expressions 
depend  on  the  unknown  n-layer  width  W  and  the  avalanche  zone  transit 
time  t^.  Solutions  for  W  and  t^  will  now  be  obtained. 

The  solution  for  W  is  based  on  the  fact  that  the  total  distance 
traversed  by  the  trailing  edge  of  the  plasma  hole  packet  from  the  formation 
of  the  plasma  to  the  start  of  residual  extraction  is 


d 

P 


“n  + V 


(72) 


d  consists  of  3  components  as  given  in  Eqs.  (36),  (42),  and  (44) 
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Using  E  =  E  „  , .  ,  we  have 

6  e  p2+bf 


E  {n  +  u  )  ?  7 

W  =  At2(2+Pbf)  130  +  8ba  -  2ba  '  5bf  '  V >' 


(74) 


Substituting  for  At  from  Eq.  (11)  we  have 
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where 
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and 


a,  = 
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(77) 


W  is  an  explicit  function  of  b^  and  E^  and  the  dependence  of  W 

on  E^  can  be  eliminated  as  follows:  From  (69)  and  using  the  fact  that 

E  «E  ,  E  may  be  written  as 
e  mm 
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Using  (63) 
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where 
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where 
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Substituting  (82)  into  (75)  gives 
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Now  return  to  the  solution  of  b^.  By  definition 
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Now 
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where 
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Using  Eq.  (30)  we  have 
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The  result  is  a  cubic  equation  in  which  can  be  easily  solved  for  its  real, 
positive  root: 


bf3(l  +k2)  +  bf2(4  +  7k2)  +  bf(4  -  2k  -  [2C  -  10]k2) 
-(4k1  +  4Ck2)  =  0 


(94) 


This  completes  the  TRAPATT  analysis. 
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J.  Numerical  Results 

This  section  gives  the  results  of  a  sample  calculation  using  the 
analysis  described  above.  The  diode  parameters  are  assumed  identical  to 
those  listed  in  Table  I  with  the  following  exceptions: 

a.  The  normalized  average  current  density  during  the  ava¬ 
lanche  is  j  =  2.  2. 

b.  The  n-layer  width  W  is  not  specified  as  an  input  parameter. 

c.  The  peak  current  Jq  is  not  specified  as  an  input  parameter. 

d.  The  normalized  time  from  the  start  of  the  current  fall  to 

the  start  of  residual  extraction  is  b  =  0.4. 

a 

The  results  of  the  calculation  are  given  in  Table  III.  The  computed 
voltage  waveform  for  one  cycle  is  given  in  Fig.  9. 

It  is  seen  that  the  junction  field  does  not  exceed  the  static  breakdown 

value  of  3  x  10  v/cm  except  during  avalanche  formation.  The  n-layer  width 

is  7.  496/im,  the  frequency  is  1.  49  GHz,  and  the  peak  current  density  is 
3  2 

5.  33  x  10  A/cm  .  These  results  may  be  compared  with  the  results  obtained 

by  applying  the  C.I.  N.  square  wave  analysis  to  a  diode  of  the  same  width 

and  assuming  the  same  peak  current.  The  outcome  of  that  calculation  is 

summarized  in  Table  IV.  The  resulting  frequency  is  948  MHz.  Because  of 

3  2 

the  large  peak  current  density  of  5.  33  x  10  A/ cm  ,  premature  avalanche 
occurs  in  this  diode  during  plasma  extraction  despite  the  assumption  of  an 
ideal  square  wave  current  drive. 

Harmonic  impedances  are  obtained  by  Fourier  analyzing  the  voltage 

waveform%of  Fig.  9.  Table  V  gives  the  calculated  impedances  through  the 

-3  2 

7th  assuming  a  junction  area  of  10  cm  .  The  efficiency  and  average  voltage 
are  also  given.  In  the  2nd  column  of  Table  V  these  values  are  compared  with 
the  corresponding  results  obtained  by  analysis  of  the  voltage  waveform  for 
the  square  wave  case  of  Table  IV.  It  is  seen  that  the  7th  harmonic  impedance 
is  positive  in  the  square  wave  case.  This  means  that  an  active  termination 
would  be  required  for  an  exact  match  at  this  frequency.  However,  as  we 
have  seen  in  Fig.  4,  the  elimination  of  the  7th  harmonic  current  in  the  tr ape- 
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DIODE  VOLTAGE  (volts) 


1/T 

Fig.  9.  One  cycle  of  voltage  waveform  computed  for 
TRAP  ATT  oscillator,  driven  by  trapezoidal  current. 
Diode  and  operating  parameters  for  this  case  are 
given  in  Table  III.  Current  rise  time  is  0.  2T. 
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TABLE  III 


A 


Summary  of  Trapezoidal  Current  Calculation  for  TRAPATT  Oscillator 


Diode  Parameters 

Same  as  in  TABLE  I  except  N-Layer  Width  (W):  Not  Specified 

Input  Parameters 

Normalized  Current  Density  During  Avalanche  (j  ):  2.  2 

cl 

Normalized  Time  to  Residual  Extraction  (b  ):  0.4 

SL 

Peak  Current  Density  (Jq):  Not  Specified 


Junction  Field  Calculations 


Start  of  Plasma  Formation  (E  ):  3.  29  x  10  v/ cm 

c  5 

Maximum  During  Plasma  Formation  (E  ):  3.  83  x  10  v/cm 

pk 

End  of  Plasma  Formation  (E  ):  1.  20  x  lO^v/cm 

p 

End  of  Current  Rise  (E  ):  1.  82  x  lO^v/cm 

e 

Start  of  Current  Fall  (Ef  ):  2.65  x  10^v/cm 

fP  5 

Maximum  During  Residual  Extraction  (E  ):  2.83  x  10  v/cm 
6  max 

End  of  Residual  Extraction  (En):  2.47  x  105v/cm 
Holding  Field  (Em):  2.99  x  105v/cm 


Other  Results 
N-Layer  Width  (W):  7.496  x  10~^cm 

3  2 

Peak  Current  Density  (Jq):  5.  33  x  10  A/ cm 
Normalized  Time  for  Plasma  Formation  (b^):  0.642 
Frequency  (f):  1.49  GHz 
Average  Voltage:  115.  6v 
Efficiency:  48.  8% 
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TABLE  IV 


Summary  of  Square  Wave  Current  Calculations  for  ». 

TRAP  ATT  Oscillator  of  Table  III 

> 

Diode  Parameters 

Same  as  Table  I  except:  W  =  7.496  x  10  cm 

Input  Parameters 

J  =  5.  33  x  103A/cm2 
o 

j  =  3.  331 
Ja 

Junction  Field  Calculations 

E  =  3.  63  x  103v/cm 
c 

E  .  -  4.  31  x  10^v/cm 

pk 

3 

E  =  1.04  x  10  v/cm 
P 

Frequency 
f  =  0.  948  GHz 
Average  Voltage 

129.  0v 
Efficiency 

52.  9%  j 
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TABLE  V 


Impedances  of  TRAPATT  Oscillators  of  Tables  III  and  IV 

-3  2 

j  Assuming  Junction  Area  of  10  '  cm 


Harmonic 

Number 

Impedance  (Ohms) 
(Trapezoidal) 

Impedance  (Ohms) 
(Square  Wave) 

1 

-30. 1  + j  0.  286 

-33.  6  -  j  2.  98 

3 

-37.  4  +  j  32.  8 

-21.  7  +  j  8.  9 

5 

0.  0 

-12.  0  +  j  14.4 

7 

o 

O' 

o 

CM 

4- 

-6. 7  +  j  16.  32 

zoidal  waveform  has  little  effect  on  the  shape  or  rise  time  of  the  waveform. 

Thus,  the  7th  harmonic  may  be  terminated  in  an  open  circuit  without  affecting  the 
trapezoidal  analysis  in  a  first-order  way,  leaving  a  requirement  for  complex 
impedance  terminations  only  at  the  fundamental  and  third  harmonics.  All 
other  harmonics,  both  even  and  odd  are  terminated  in  open  circuits.  From 
Table  V  it  is  seen  that  the  fundamental  impedances  in  both  cases  are  approxi¬ 
mately  pure  real  and  differ  in  magnitude  by  only  about  10%.  In  the  trape¬ 
zoidal  case  the  real  part  of  the  3rd  harmonic  impedance  is  nearly  double 
the  corresponding  impedance  in  the  square-wave  case.  The  reactive  component 
is  more  than  a  factor  of  three  larger.  The  larger  3rd  harmonic  impedance- 
results  principally  because  the  3rd  harmonic  of  the  trapezoidal  current  wave¬ 
form  is  twice  that  of  the  square  wave. 

The  magnitude  of  the  3rd  harmonic  impedance  is  also  higher  than  the 
magnitude  of  the  fundamental  impedance  in  the  trapezoidal  case,  making  it 
easier  to  realize  the  3rd  harmonic  termination  experimentally. 

The  efficiency  in  the  square-wave  case  is  52.  9%  as  compared  to  48.  8% 
in  the  trapezoidal  case.  This  approximate  6%  difference  in  efficiency  occurs 
primarily  because  the  ratio  of  fundamental  to  dc  in  the  trapezoidal  current 
waveform  is  approximately  6%  less  than  the  same  ratio  for  the  square  cur- 
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rent  waveform,  while  the  resulting  ratios  of  fundamental  to  dc  voltage  are 
identical  to  within  1%  in  the  two  analyses. 

K.  Conclusions 

An  approximate  analysis  of  TRAPATT  operation  has  been  described 
which  employs  physically  realizable  waveforms  and  which  results  in  realistic 
conditions  on  the  oscillator  harmonic  terminations.  This  analysis  is  based 
on  specific  assumptions  concerning  the  shape  of  the  driving  current  wave¬ 
form  and  the  timing  of  the  critical  transitions  in  the  diode.  It  is  thus  more 
of  an  existence  demonstration  than  a  general  analysis  of  TRAPATT  operation. 

The  results  demonstrate  the  realizability  of  continuous  high  efficiency 
oscillations  in  an  idealized  avalanche  diode  structure  using  a  trapezoidal 
current  waveform.  The  calculations  also  indicate  that  while  a  previous 
TRAPATT  analysis  assuming  an  ideal  square  wave  current  gives  a  reasonable 
approximation  to  the  fundamental  impedance  and  efficiency,  the  3rd  and  5th 
harmonic  impedances  are  significantly  altered  in  going  from  the  unrealizable 
square  wave  to  the  realizable  trapezoidal  current  drive. 

One  of  the  principal  reasons  for  developing  a  finite  rise  time  analysis 
was  to  better  understand  and  prevent  premature  avalanche  at  the  end  of  the 
plasma  extraction  period.  The  results  of  this  trapezoidal  analysis  indicate 
that  premature  avalanche  can  be  avoided  when  using  a  finite  fall  time  current 
drive  if  the  current  fall  begins  early  in  the  plasma  extraction  process.  How¬ 
ever,  it  is  clear  that  precise  waveform  control  is  required  to  accomplish 
this,  and  that  small  variations  in  operating  frequency  or  harmonic  terminations 
could  result  in  major  changes  in  the  operating  mode  of  the  oscillator.  Thus, 
an  oscillator  operating  in  this  finite  rise-time  or  trapezoidal  current  mode 
will  be  inherently  narrowband  and  will  operate  only  over  a  small  range  of 
3rd  and  5th  (in  this  case  the  5th  is  an  open  circuit)  harmonic  impedances. 
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